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Abstract 


This  report  discusses  some  further  developments  of  an  analytical 
solution  to  the  3-D  Griffith  crack  problem.  The  analysis  shows  the 
stresses  at  the  corner  points  to  be  singular  of  the  order  (1/2  +  2v). 
Moreover,  the  stress  boundary  conditions  at  the  plate  faces  are  shown  to 
be  proportional  to  (h  -  z),  at  the  upper  face,  and  to  (h  +  z),  at  the 
lower  face. 


debility  in  current  fracture  mechanics  work  is  the  ignorance  of 
the  effects  of  thickness  on  the  mechanism  of  failure.  For  example,  the 
common  experimental  observation  of  a  change  from  ductile  failure  at  the  edge 
to  brittle  fracture  at  the  center  of  a  broken  sheet  material  has  so  far 
defied  analysis.  Moreover,  fracture  toughness  of  a  material  is  presently 
assumed  to  be  constant.^Yet  engineering  tables  prescribe  two  different  ^ 
values  for  thin  and  thick  plates.  More  frustrating  yet  is  the  fact  that 
there  exists  no  precise  definition  as  to  what  represents  a  'thin'  or  a 
'thick'  plate.  Obviously,  the  concept  of  fracture  toughness,  as  it  is 
presently  defined,  is  an  inappropriate  parameter  to  use  for  the  prediction 
of  failures  due  to  fracture.  Yet  an  orderly  theoretical  attack  on  the 
problem  can  provide  important  guidance  to  these  and  other  phases  of 
fracture  research.  The  most  potent  mathematical  tool  for  this  attack  is  the 
linear  theory  of  infinitesimal  elasticity  as  applied  to  a  cracked  plate  of 
finite  thickness.  Although  this  theory  cannot  include  the  nonelastic 
behavior  of  the  material  at  the  crack  tip  per  se,  it  can  evince  many 
characteristics  of  the  actual  behavior  of  a  cracked  plate,  including  those 
due  to  thickness. 

For  example,  such  information  could  be  valuable  to  the  understanding 
and  the  solution  of  the  corresponding  3-D  elastoplastic  problem.  While  it  is 
true  that  this  theory  cannot  give  us  the  exact  stress  state  at  a  point  in  the 
interior  of  the  shear  lip,  it  can,  however,  prescribe  fairly  accurately  the 

shear  lip  envelope.  Thus  the  theory  of  linear  elasticity  is  a  logical 

★ 

fountainhead  for  detailed  theoretical  study. 

★ 

It  is  for  these  reasons  that  a  National  Workshop  on  Three-Dimensional 
Fracture  Analysis,  held  at  Battelle  on  April  26  -  28,  1976,  identified  this 
as  one  of  the  Benchmark  problems  in  the  field  of  Fracture  Mechanics. 


Introduction 


A  major 


The  mathematical  difficulties,  however,  posed  by  three-dimensional  crack 


problems  are  substantially  greater  than  those  associated  with  plane  stress 
or  plane  strain.  As  a  result,  there  exist  in  the  literature  very  few 

analytical  papers  that  deal  specifically  with  the  three-dimensional  stress 

★ 

character  at  the  base  of  a  stationary  crack.  Simultaneously,  in  the  last 
two  decades,  numerous  attempts  have  been  made  to  obtain  a  finite  element 

solution  by  very  capable  researchers.  Unfortunately,  they  too  experienced 

★  ★ 

difficulties  and  their  respective  results  were  contradictory. 

Review  of  author' s  past  work . 

In  1973,  Folias,  using  a  method  developed  by  the  Russian  elastician 
Lur'e  [1]  and  the  application  of  Fourier  Integral  Transforms,  constructed 
a  solution  [2]  to  Navier's  equations  for  a  mixed  boundary  value  problem, 
that  of  a  3-D  Griffith  crack  (see  fig.  1).  The  integral  representations 
were  subsequently  expanded  asymptotically  in  the  inner  layers  of  the  plate 
and  the  displacement  and  stress  fields  were  found  to  be: 

(i)  Displacements: 
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For  a  historical  discussion  see  reference  [3]. 

**Refer  to  proceedings  of  Workshop  on  Three-Dimensional  Fracture  Analysis 
held  at  Battel le,  1976. 
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where  /.  is  a  function  of  Poisson's  ratio  v  and  c/h. 

The  above  solution  reveals  the  following  important  characteristics, 
which  are  applicable  orly  in  the  inner  layers  of  the  plate: 

(1)  the  stresses  possess  the  usual  1//T  singularity, 

(2)  the  stresses  possess  the  usual  angular  distribution, 

(3)  the  stress  intensity  factor  Kj  is  a  function  of  z, 

(4)  exact  plane  strain  conditions  exist  only  on  the  plane  z  =  0, 

★ 

(5)  a  pseudo  plane  strain  state  exists  and  the  equation 

- Z  =  V('X  +  cy)  (10) 

(6)  as  the  plate  thickness  2h  — -  ®,  the  plane  strain  solution  is 
recovered, 

(7)  as  Poisson's  ratio  v  —  0  ,  the  plane  stress  solution  is 
recovered. 


The  reader  should  note  that  the  stresses  are  functions  of  z,  therefore 
the  term  'pseudo'  plane  strain. 


and  also  serve  as  check  points  for  the  validity  of  the  solution,  at  least 
as  v  — *■  0  and  as  h  — *  ®. 

One  of  the  great  advantages  of  an  analytical  expression  is  that  it  can 
easily  be  verified  as  to  whether  it  represents  a  solution  to  a  problem  by 
direct  substitution  into  the  governing  equations  and  the  appropriate 
boundary  conditions.  The  reader  should  note  that  the  above  displacement 
field  satisfies  all  Navier's  equations  and  that  the  stress  field  satisfies 
the  following  boundary  conditions 

Txy  =  Tyz  =  0  ;  for  *  =  0  and  r 


for  $ 


Inasmuch  as  the  solution  represents  an  asymptotic  expansion  which  is 
valid  only  in  the  inner  layers  of  the  plate,  naturally  the  stress  field  can 

not  be  expected  to  satisfy  the  boundary  conditions  on  the  plate  faces,  i.e. 

★ 

at  2  =  x  h  .  The  author,  however,  shows  in  reference  [2]  that  in  the 

neighborhood  of  the  corner  point  additional  terms  also  contribute  to  the 

★★★ 

same  order  of  singularity  and  therefore  must  be  accounted  for.  Moreover 
it  was  shown  that  all  stresses  there,  inclusive  of  and  x^,  possess  the 
same  order  of  singularity.  He,  subsequently,  ventured  to  examine  the 
character  of  the  stress  singularity  which  prevails  at  the  corner  point,  i.e 
the  point  where  the  crack  front  meets  the  free  surface  of  the  plate.  After 


See  pages  669  -  671 . 

This  matter  will  also  be  addressed  at  a  later  time. 

It  should  be  pointed  out  that  the  integral  representations  do  satisfy  the 
boundary  conditions  there  and  that  this  comment  refers  only  to  the  recovery 
of  the  complete  local  expansion  from  the  integral  representations. 


considerable  effort,  it  was  concluded  that  the  displacements  are  proportional 

to  and  similarly  the  stresses  are  proportional  to  p-l/2-2v. 

Researchers  in  the  field  of  fracture  mechanics,  however,  were  reluctant 

★ 

to  accept  at  that  time  the  possibility  of  an  infinite  displacement  field  , 
on  the  basis  of  physical  intuition.  Consequently,  the  results  were  con¬ 
sidered  to  be  controversial  and  the  following  legitimate  and  fundamental 
questions  were  raised: 

1.  Is  the  solution  of  this  notorious  difficult  problem  unique  and  if  so 
under  what  conditions? 

2.  Is  an  infinite  displacement  field,  in  this  case,  admissible? 

r  21 

3.  Does  the  symbolic  method  adopted  by  the  author1-  J  generate  a 
'complete'  set  of  eigenfunctions  for  the  solution  of  Navier's 
equations? 

★  ★ 

The  answers  to  the  above  questions  were  given  by  Prof.  Calvin  Wilcox. 

First  of  all,  he  was  successful  in  proving  [4]  that  a  displacement  field 
which  satisfies  the  condition  of  local  finite  energy  is  unique.  This  of 
course  is  quite  a  departure  from  our  traditional  2-D  fracture  mechanics 
thinking,  for  the  displacements  may  now  be  allowed  to  be  singular.  Conse¬ 
quently,  one  may  not  a  priori  assume  them  to  be  finite  as  it  is  customarily 
done.  In  general,  such  an  assumption  makes  the  class  of  solutions  too 


Infinite  displacement  fields  are  not  uncommon  in  the  theory  of  linear 
elasticity.  For  example,  it  is  well  established  that  the  Boussinesq's, 
Kelvin's  and  Cerruti's  problems  have  displacement  fields  that  behave 
like  ui  -  The  reader  should  also  be  reminded  that  the  reported 

results  of  the  crack  problem  correspond  to  the  complementary  solution, 
i.e.  the  interior  of  the  crack  faces  is  subjected  to  a  uniform  fluid 
pressure. 

This  work  was  carried  out  at  the  request  of  the  author  and  was  supported 
by  AFOSR  Contract  F4. 962077. C0053  which  the  author  had  at  the  time. 


restrictive  and,  as  a  result,  one  may  not  find  the  complete  solution  to  the 
problem.  On  the  other  hand,  the  solution  could  very  well  give  finite  dis¬ 
placements  everywhere!  Be  that  as  it  may,  physical  intuition  should  be  used 
with  extreme  caution. 

★ 

Second,  he  was  able  to  show  [5]  that  the  Fourier  integral  expressions 
representing  the  general  solution  to  Navier's  equations  are  complete  and, 
furthermore,  the  'symbolic  method'  used  is  justifiable.  In  order  to  prove 
this,  he  used  a  double  Fourier  integral  transform  in  x  and  y  and  subsequently 
a  contour  integration  to  recover  precisely  the  same  expressions  as  those 
reported  by  Fo 1 i as  in  [2]. 

Finally,  it  remains  to  determine  explicitly  the  stress  field  ahead  of 
the  crack  tip  and  throughout  the  thickness  of  the  plate.  Moreover,  it 
would  be  interesting  to  see  if  the  asymptotic  expansion  in  the  inner  layers 
of  the  plate  matches  the  stress  field  given  by  equations  (4)-  (9). 


The  Complementary  Displacement  and 
Stress  Fields 

From  reference  [3],  we  have  the  following  complementary  displacement 
and  stress  fields. 


( 1 )  the  displacement  fie l_d_: 
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See  Eqns.  (52)  -  (54)  of  [4]. 
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where  for  simplicity  we  have  adopted  the  following  definitions: 
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with  A  and  B  as  functions  of  t  and  an  respectively. 


the  roots  of  the  transidental  equation 
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the  stress  field: 
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By  direct  substitution,  it  can  easily  be  ascertained  that  the  fore¬ 
going  complementary  displacement  field  satisfies  Navier's  equations  and 
furthermore  the  corresponding  stresses  o7^C\  r  t  do  vanish  at 

l.  tA  y  c. 

the  plate  faces  z  =  r  h.  The  remaining  boundary  conditions: 
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‘  0 


m  +  1  4  ,..2  . ,  }y=0  "  ‘  ^  2G  *  ;  'x>  <  c  * 


s*  ? |y  i 


Next,  upon  utilizing  eq.  (30)  into  eq.  (31),  one  finds:  after  some 
rearrangement: 


3H 


~2  ‘  ~2  *  lT7)y=0[2("!'1)cos(!vh,cos(evz) 


m  i  h  s  i  n  ( £  h )  cos  ( f  z  )  +  m  f  z  cos  ( r  h )  sin  ( r  z )  ] 

V  V  V  V  V  V 


8m 


m-  2  ^  Av  (T7Jy=0  J 


( 


-  2  “  3H 

i,  Av(l7>y0[r'vhSln(^h)CO5lV) 


oz 


B  z  cos  (r,  h)  sin  (B  z)  ]  } 

v  J  V 


3 

1  3  I  3 1  o  O  o  3  I 

.  {  .  2  (  m~  1  )  —1  -  2  —  +  2  7 2  _ l  } 

'  £  '  m  4 .  1  /  Z  3  v  -  “  1  Z  ^  J 


m+  1  '  3x  3x  m  +  1  L  3  '  y=0  ’ 

3x  3 


for  all  x  . 


The  problem,  therefore,  has  been  reduced  to  that  of  solving  equations 


(30),  (32)  and  (33)  for  the  complex  coefficients  A  ,  the  real  coefficients 

V 

B  ,  and  the  respective  functions  H  (x,y;6  )  and  H  (x,y;a  )  subject  to  the 
n  v  v  n  n 

continuity  conditions 


?H  , 

^  I  n 
&y  ;y=0 


^y  !:y"0 


X  >  c 


X  !  >  c 


(34) 


(35) 


However,  before  one  engages  into  a  lengthy  numerical  procedure  for  the 

solution  of  the  system:,  it  would  be  desirable  to  construct,  if  possible,  the 

★ 

functions  H  (x,y;r  )  and  Hn(x,y;r,  )  analytically.  After  considerable 
mathematical  manipulations  it  is  possible  to  construct  these  functions. 

it 

For  example,  without  going  into  the  mathematical  details  ,  H*  is  found  to 
be : 


=  -  {  x  e 


"Vi 


2  ✓ 


[/r,  +  |  y  |  +  /r. 


TTQ. 


+  — ^zz-iy r2 +  iy ! 

2  *  in 


/r. 


y  I  ]  e 


anr2 


+  4(x-p-^){e  n  Erf[/a(r1-|yi)]  (36 

L  ^an  n  1 

+  e  Vy'  Erf  C/a  (rn  +  |  y  j)  ]  +  + 

“n 

[  e  n  'y'  Erf  [ /un(r2  -  |y|)]  -e^  '  *  '  Erf  [/  ^(i^  +  |y \)V 


*See  footnote  on  next  page. 


where  we  have  adopted  the  following  definitions 


A  similar  expression  may  be  obtained  for  H^(x,y;Bv)  as  a  combination  of 

the  function  (36)  and  its  derivatives.  This  step  was  very  important  for  it 

enables  one  to  extract  the  dependence  in  x  and  y  analytically  rather  than 

numerically.  In  summary,  therefore,  the  problem  has  now  been  reduced  to 

that  of  solving  equations  (30),  (32)  and  (33)  for  the  complex  constant 

coefficients  A  and  the  real  but  constant  coefficients  Bn-  Moreover,  the 

displacement  and  stress  fields  may  now  be  expressed  in  terms  of  one  summation 

rather  than  one  summation  and  one  integral. 

Although  the  numerical  solution  of  these  coefficients  is  now  possible, 

from  which  the  singularity  may  then  be  determined  (at  least  in  principle), 

we  believe  that  such  procedure  may  not  lead  us  to  all  the  singularities  of 
★ 

the  problem  and  that  such  information  should  be  extracted  by  analytical 
means . 

The  behavior  of  the  solution  at  the  corner . 

One  of  the  acadamic  questions  that  has  defied  researchers  for  a  long 
time  is  the  order  of  the  singularity  that  prevails  at  the  corner  point. 

Thus,  in  order  for  us  to  study  the  stress  distribution  in  the  vicinity  of 
that  point,  we  introduce  the  following  function 

f  h  (h  +  z )  =  I  A  (  )y=0  cos  [  bv  (h  +  z)  ]  ,  (39) 

in  view  of  which  the  homogeneous  portion  of  equation  (33)  may  now  be 
written  as 


★ 


This  is  because  numerical  procedures  tend  to  average  local  affects 
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where 


E 


E 


-  6  (1  +  3  vQ) 


12vq  +  20  vQ  +  21/4 


-2v0  (l  +  6v0) 


-  6  (1+  v0) 


4  (2+3  .0) 


m  =  -2(1  +  3.0)  . 


(45)  -  (51) 


Finally,  the  solution  of  equation  (44)  can  be  found  by  the  application  of 
the  method  of  Frobenious,  i.e. 


F  ( ' )  =  cr 

n=0 


Without  going  into  the  details,  the  indicia!  equation  is  found  to  be 


(r  -  1)  (  r  -  2)  +  (6  ,n+2)  S(f  -  1)  ♦  2vn  (1  +  6  vn)  S 


+  4v0M2v0-l)  =  0  , 


the  roots  of  which  are: 


T  -2  double  root 


E 


One  concludes,  therefore,  that 


f(h)(h+  2) 


Cn  (c  -  x) 


•1/2 


{  (c  -  x)  +  (h -  z ) ‘ 


1-v 

)  H(n) 


(55) 


where  H ( r> )  represents  an  even  power  series  in  n  the  behavior  of  which  for 

small  values  of  n,  i.e.  around  the  point  ri  =  0,  is  shown  in  fig.  2.  A 

★  ★ 

similar  asymptotic  expansion  around  the  points  n  =  +  1  shows  that 


H( • )  =  (l—,2)’  I  CM-  r2)n 

n=0  n 


(56) 


3 

with  ■,  =  -^  .  The  behavior  of  H  around  the  point  n  =  -l,  for  example,  is 

shown  in  fig.  3.  It  should  be  emphasized  that  there  exist  no  singular 

points  of  any  kind  between  zero  and/one.  Hence  the  function  H(n),  in 

between,  represents  a  smooth  curve.  Pending  a  further  analytical  or 

numerical  study  of  H(',)  the  function  has  been  approximated  by  the  dotted 
★  ★★ 

curve. 

In  view  of  eq.  (55),  the  function  of  f^(h+z)  has  now  been 
determined.  While  there  are  many  ways  that  one  can  now  use  to  complete  the 
rest  of  the  problem,  the  authors  suggest  the  following  method  in  order  to 
extract  an  expli ci t  analytical  form  of  the  solution.  First,  one  needs  to 
construct  a  harmonic  function  v(x,y,h+z)  such  that 


★ 

The  reader  should  note  that  this  is  only  one  out  of  the  three 
solutions  and  that  the  double  root  leads  to  a  logarithmic  solution. 

The  reader  should  note  that  the  other  two  values  of  y  =  0  and  >  =  ^ 
correspond  to  the  logarithmic  and  the  6=1-  2vq  solutions, 
respectively. 

In  the  event  that  a  numerical  scheme  is  used,  the  matching  of  the  two 
curves  must  be  done  with  care. 


Discussion  of  the  results 


In  order  for  us  to  see  the  form  of  the  singularity  which  prevails  at 
the  vicinity  of  the  corner  point,  we  introduce  the  local  spherical 
coordinates  defined  by  the  transformation: 
c  -  x  =  p  sin  r  cos  c 

y  =  e  sin  e  sin  c  (59) 

h  -  2  =  p  cos  6 , 

which  lead  to* 

2 

f(h)  (h  +  z)  =  CQr  2  '  2v  — ! -  H  (cos  e)  ,  (60) 

/sin  b 


frorr  which  it  is  evident  that  the  displacements  are  proportional  to 

1/22  1/2  2  v 

c  7  and  the  stresses  proportional  to  c"  ~  v.  We  observe,  there¬ 

fore,  that  this  is  the  same  result  as  that  reported  in  reference  [2]. 

Let  us  next  examine  the  behavior  of  the  function  f  in  the  inner  layers 
of  the  plate,  i.e.  for  all  x  -  c  =  €  <<  h  -  2 


:(h)  ,u  ,  .  "0 


(h  +  2)  =  -Mh 

a 


z) 


2-2v 


{  1  + 


( h  -  2 ) 1 


1-v 

)  H  ( 


(h  -  2)- 


(61) 


=  ^  (h  -  z)2*2v  H  (1)  ♦  0  ( €° )  , 

Jl 

which  again  leads  to  the  same  type  of  results  as  those  reported  in  ref.  [2], 
We  now  have,  therefore,  a  better  understanding  of  the  type  of  expansions 
that  equations  (4) -  (9)  represent. 


* 

The  reader  should  recall  that  the  function  f 


(b) 


is  defined  on  the  plane 


It  would  be  appropriate  at  this  time  to  reexamine  the  question  as  to 


whether  boundary  conditions  and  can  indeed  be  satisfied  at  the 

★ 

plate  faces.  An  examination  of  these  stresses  shows  that  they  can  be 


expressed  in  terms  of 
„2 

(h  -  z)  H 


a"d  fz 


and  similar  expressions.  The  former  is  automatically  satisfied  by  virtue 

of  the  factor  (h-z).  The  latter,  after  the  odd  derivative  is  taken,  a 

factor  (h-z)  is  kicked  out  thus  forcing  it  to  vanish.  The  reader  should 

recall  that  f^  is  an  even  function  of  (h-z).  All  other  terms,  which  are 

(P) 

of  a  lesser  order,  are  negated  by  the  particular  solution  [fv  One 

concludes,  therefore,  that  such  an  expansion  can  indeed  satisfy  the 
boundary  conditions  at  the  plate  faces  and,  as  a  result,  is  admissible. 

It  would  be  interesting  next  to  compute  the  displacement  function 
v(x,o,z)  in  the  neighborhood  of  the  corner  point.  Without  going  into  the 
mathe-.atical  details  one  finds  after  some  rearrangement  that 


v<c>  =  CQ  .  c  -  x  t'2'  G(n)  + 


(62) 


where 


G(-)  =  [(A  VQ2  -  4Vq)  n2  +  (-  4  vq2  +  vq)]  H  + 


+  [(3  -  4  v 0 ) r, 3  +  (4vq  -  |)n]H’  +  (n2  -  D  n2  H". 


(63) 


*The  reader  should  note  that  there  was  never  a  question  as  to  whether  the 
solution  in  the  integral  representation  form  did  just  that. 

**  /  \ 

The  reader  should  also  take  into  account  eq.  (30). 

The  ...  here  indicate  terms  of  higher  order. 


A  plot  of  the  function  G(s)  is  given  in  Fig.  3  for  various  Poisson's 
ratios.  The  value  of  n  =  0  corresponds  to  the  free  surface  while  n  =  -1 


corresponds  to  the  crack  front.  The  constant  Cq  is  proportional  to  the 
applied  stress  and  a  function  of  Poisson's  ratio.  Thus,  although  the 
amplitude  of  v  may  change,  the  trend  will  remain  the  same.  An  examination 
of  the  curve  also  shows  that  as  one  moves  from  the  crack  front  to  the  free 
surfaces  of  the  plate  the  crack  widens,  a  result  that  meets  our  expecta¬ 
tions. 

Notice  also  that  for  zero  Poisson's  ratio  the  function  G  is  constant 
and  the  displacement  now  becomes 

lin  v(c)  (* ,o , z )  =  {  lim  Cn  •  const.  /  c  -  x  .  (64) 

v-*0  v-C 


This  result  also  meets  our  expectations  for  it  represents  an  exact 
solution. 

Let  us  next  examine  the  behavior  of  the  displacement  in  the  neighbor¬ 
hood  of  the  cracV  front,  i  .e.  when  r  =  -1.  The  function  H(r,)  then  behaves 
like 

H ( •  )  -  (1  -  -2)3/4  , 
which  then  leads  to 


v(c»-.  const.  .W'.TO-||.(VW. 


cons 


t.  .w-2  .i «  ••• 


(65) 


See  eq.  (56). 
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Conclusions 


-  i 

As  a  result  of  this  study,  the  authors  draw  the  following  conclusions. 

First  of  all  the  displacements  in  the  neighborhood  of  the  corner  point 
are  at  least  proportional  to 

ui  x  ;1/2'2v  fi (e ,c)  (66) 

and  that  the  functions  f ^ (6,4)  do  not  vanish. 

Second,  while  it  is  true  that  for  Poisson's  ratios  greater*  than  1/4 

the  displacements  are  infinite,  the  local  strain  energy  however  remains 

finite.  Consequently,  such  a  solution  is  admissible  and,  if  it  represents 

r  41 

'a  solution',  it  is  the  only  solution1-  J.  Infinite  displacements  are  not 
new  to  the  field  of  linear  elasticity.  For  example,  Cerruti's  problem 
leads  to  higher  order  singularities  (displacements  <v  l/c).  It  should  be 
emphasized  that  the  problem  which  the  authors  have  solved  is  that  of  the 
complementary  3-D  Griffith  crack,  where  the  crack  faces  now  are  subjected 
to  a  uniform  internal  pressure  Og. 

Third,  the  solution  is  not  separable  in  spherical  coordinates.  Thus 
solutions  that  do  invoke  such  an  assumption  may  or  may  not  lead  to  the 
correct  or  complete  solution  of  the  problem. 

Fourth,  physically  what  the  solution  really  shows  is  that  linear 
theory  is  inadequate  in  predicting  the  actual  behavior  of  the  material  at 
such  corner  points.  This,  however,  should  not  be  interpreted  as  if  no 
practical  information  can  be  extracted  from  such  a  solution.  An  example  of 
this  is  the  2-D  solution  which  predicts  infinite  stresses  at  the  base  of  the 
crack. 

Fifth,  the  material  cannot  sustain  such  large  stresses  at  the  corner 
and,  as  a  result,  it  must  relax.  This  relaxation  is  accomplished  in  two 

♦For  v  *  1/4,  the  displacements  are  'v  pV2-2(l/4)  .  p° 


In  this  case  u  Inp. 
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ways.  Initially  the  material  yields,  thus  forming  the  shear  lips,  and 
simultaneously  it  fractures,  thus  forming  the  well  known  'tunneling  effect'. 
The  authors  believe  that  this  fracture  initiates  first  at  the  corner, 
where  the  material  tries  to  smooth  out  the  90°  angle  and  subsequently  its 
crack  advances  in  the  manner  shown  on  fig.  4. 

Sixth,  present  experimental  results  may  not  be  used  to  prove  or  dis¬ 
prove  a  solution  for  the  following  reasons.  In  experimental  work  a  crack 
is  usually  introduced  by  mechanical  means  and  is  subsequently  fatigued  in 
order  to  attain  a  sharp  crack  tip.  The  process  of  fatigue,  however,  will 
undoubtedly  smooth  out  these  sharp  corners  thus  altering  the  geometry  of 
the  problem.  Moreover,  if  the  displacements  are  indeed  very  large,  the 
material  will  yield  and  the  shear  lips  will  be  formed.  In  such  a  case, 
the  principle  of  superposition  is  no  longer  valid  and  any  such  comparison 
may  be  meaningless,  unless  the  load  in  the  experimental  procedure  is 
applied  to  the  crack  faces  in  which  case  one  then  compares  'complementary' 
problems . 

Seventh,  it  has  been  observed  experimentally  that  on  the  free  surfaces 
of  a  plate  and  around  the  neighborhood  of  the  crack  tip,  a  small  dip 
appears  during  the  process  of  deformation.  This  dip  substantiates  the  fact 
that,  whenever  a  load  is  applied  to  the  plate,  it  forces  the  material  to 
smooth  out  the  ninety  degree  corner.  Moreover,  the  larger  the  load  is  the 
more  noticeable  the  dip  is.  As  in  the  2-D  case,  large  forces  of  atomic  or 
molecular  attraction  of  the  order  of  the  "theoretical  strength"  prevail  at 
such  corners.  These  "cohesive  forces",  attempting  to  smooth  out  the  right 
angle,  pull  the  material  towards  the  center  of  the  plate  thus  forcing  the 
crack  front  to  'buckle'  and  tunnel  its  way  into  the  center  of  the  plate. 


Eighth,  the  authors  believe  that  there  exist  two  different  types  of 


solution  expansions  that  intermix  to  form  the  various  singularities  of 
the  problem.  This  phenomenon,  coupled  with  the  fact  that  the  solution 
is  not  separable  in  spherical  coordinates  may  offer  an  explanation  as  to 
why  other  authors  obtained,  numerically,  results  of  a  lesser  order 
singularity. 

Engineering  fracture  mechanics  can  deliver  the  methodology  to  com¬ 
pensate  the  inadequacies  of  conventional  design  concepts  that  were  based 
on  tensile  strength  or  yield  strength.  While  such  criteria  are  adequate 
for  many  engineering  structures,  they  are  insufficient  when  the  likelihood 
of  cracks  exists.  Now,  after  approximately  three  decades  of  development, 
fracture  mechanics  has  become  a  useful  tool  in  design,  particularly  with 
high  strength  materials.  Interestingly  enough,  this  advancement  has  been 
based  primarily  on  the  2-D  solution  of  the  Griffith  crack  problem.  We 
believe  that  the  3-D  solution  of  the  Griffith  crack  problem  will  not  only 
enable  us  to  understand  the  mechanism  of  fracture  propagation  better,  but 
it  will  also  expand  our  horizons  for  future  research  and,  hopefully,  will 
contribute  to  the  advancement  of  the  field  of  fracture  mechanics  to  higher 
levels  of  safe  design  against  fracture.  It  is  for  these  reasons  that  this 
study  was  undertaken. 
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